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Abstract. The present work is devoted to the study of dynamical features of 
Bohmian measures, recently introduced by the authors. We rigorously prove 
that for sufficiently smooth wave functions the corresponding Bohmian mea- 
sure furnishes a distributional solution of a nonlinear Vlasov-type equation. 
Moreover, we study the associated defect measures appearing in the classical 
limit. In one space dimension, this yields a new connection between mono- 
kinetic Wigner and Bohmian measures. In addition, we shall study the dy- 
namics of Bohmian measures associated to so-called semi-classical wave pack- 
ets. For these type of wave functions, we prove local in-measure convergence 
of a rescaled sequence of Bohmian trajectories towards the classical Hamilton- 
ian flow on phase space. Finally, we construct an example of wave functions 
whose limiting Bohmian measure is not mono-kinetic but nevertheless equals 
the associated Wigner measure. 



1. Introduction 

We consider the time-evolution of quantum mechanical wave functions ip e {t, •) G 
L 2 (R d ; C) governed by the Schrodingcr equation: 

(1.1) tedti; 5 = ~^A^ +V(x)r, i> e {t = 0, x) = #j € L 2 (R d ), 

where x e M. d , (eM, and V G L°°(R d ;R) a given bounded potential (satisfying 
some additional regularity assumptions given below). In addition, we have rescaled 
all physical parameters such that only one semi- classical parameter < e < 1 
remains. We shall from now on assume that H^olU 2 = lj which is henceforth 
propagated in time, i.e. 

(i-2) \mt)\\& = ui\\ L * = i. 

In addition, we also have conservation of energy, i.e. 

E e {t):=^- [ \V^ £ (t,x)\ 2 dx+ [ V(xW(t,x)\ 2 dx = E £ (0). 
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Throughout this work, we shall assume that 

(1.3) sup E £ (0) < +00. 

0<e<l 

In other words, we assume to have bounded initial energy, uniformly in e. In 
view of (|1.2[) , one can define out of ip £ {t,x) G C real- valued probability densities 
from which one computes expectation values of physical observables. Possibly, the 
two most important such densities are the position and the current- density, given 

by 

(1.4) p e (t,x) = \^(t,x)\ 2 , J%t,x) = elm^(t,x)Wr{t,x)). 

Already in 1926, the same year in which Schrddinger exhibited the eponymous 
equation, it has been realized by Madelung [TB] that these densities can be used to 
rewrite (jTTTj) in hydrodynamics! form. The corresponding quantum hydrodynamic 
system reads 




(1-5) < , ,r , „ fJ'*J<\ , p£w= ^ £v /A A 



2 

More precisely, it can be proved that under sufficient regularity assumptions on 
V(x) and ip e (t, x), each of the nonlinear terms arising in this system is well-defined 
in the sense of distributions, see [T^l Lemma 2.1]. 

The quantum-hydrodynamic system (| 1 . 5|) can also be seen as the starting point 
of Bohmian mechanics [H [3] . In this theory, one defines an e-dependent flow- map 

X £ :RxR d ^ R d ; x H. X £ (t, x) 

via the following differential equation 

X s (t, x) = u £ (t, X £ (t, x)), X £ (0, x) = x G R d , 

where the vector field u £ is (formally) given by 

, J £ (t,x) fVip £ (t,x) 
u £ (t,x) := ; ' ( = elm ' 



p £ {t,x) V ^(*>») 

and the initial data is assumed to be distributed according to the measure /5g (x) 
\i^q(x)\ 2 . Note that in terms of p £ and J £ , the kinetic energy reads 

Ekin :=^r f \Vif{x)\ 2 dx = \ [ dx + € - [ \V^\ 2 dx. 



2 y R d p e {x) 

From energy conservation, we therefore conclude u £ G L 2 (R d , p £ dx) , but not nec- 
essarily continuous. Moreover, it has been rigorously proved in [3] (see also [28]) 
that X £ (t, ■) is for all f G R well-defined p^ - a.e. and that 

p%t,x)=X £ #p £ (x), 

i.e. p £ (t,x) is the push-forward of the initial density p%(x) under the mapping 
Xf : x H> X £ (t,x), see the definition in Equation (|4.1j) below. 

Whereas, the above dynamics x i-> X £ (t, x) G K d is stated in physical space, 
Bohmian mechanics can also be reformulated on phase space M. d x R^, by using the 
concept of Bohmian measures, recently introduced by the authors in [T9"] : 
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Definition 1.1. Let e > 0. For ip £ E H 1 (M. d ), with associated densities p £ ,J £ , 
given by (|1.4[) . the Bohmian measure (3 £ = (3 £ [ip £ ] € M + (R d x R d ) is defined by 

((3 £ ,<p):=[ p £ (x)^(x,^-)dx, V^eCo^xlJ), 

where C (R^ x R d ) denotes the space of continuous function vanishing at infinity 
and M. + {K. d . x R^) the set of non-negative Radon measures on phase space. 

It has been shown shown in [TO] that if ij) £ {t,x) solves (jl.ip . then the correspond- 
ing Bohmian measure (3 £ (t,x,p) is the push- forward of 

(1.6) P e Wo) = Po( x >P) = P%{x)8(jp - ul{x)) : 
under the e-dependent phase space flow 

(1.7) $f : (x,p) i ^ (X e (i,x,p),P e (i,x,p)) 
induced by 

J X £ = P e , 

\ p £ = -w(x £ ) - wj(t,x e ), 

where Vg(t,x), denotes the so-called Bohm potential (see [TU1 ITT1 |2"5] ): 
(1.9) ^M)— 

More precisely, under mild regularity assumptions on the flow $f is shown to 
exists globally in time for almost all (x,p) £ M. 2d , relative to the measure f3 £ and 
is continuous in time on its maximal open domain, cf. |19| Lemma 2.5]. Note that 
the specific form of the initial data (|1.6j) implies that the phase-space flow <3?f is 
projected onto the graph of Uq, i.e. 

C £ :={(x,p)£RixR*:p = u e (x)}, 

whose time-evolution is governed by the Bohmian flow (|1.8|) . 



(1.8) 



The fact that f3 £ (t) = $f #/3q(x), is usually called equivariance of Bohmian 
measures (see [10]) and makes I3 £ {t) a natural starting point for the investigation 
of the classical limit as e — > 0+ of Bohmian mechanics. In [TO] we gave an extensive 
study (invoking Young measure theory) of the possible oscillation and concentration 
phenomena appearing in f3 £ as e — ¥ + and compared our findings to the by now 
classical theory of Wigner measures, cf. [TOJ [T7J [37] . One thereby associated to any 
wave function ip £ a Wigner function w £ [ip £ ] = w £ , defined by |29j : 

W £ (t, x,P):=j^ I r (t, X - £ -y)r (t, * + |y) e- dy. 

It is well known that although w £ (t,x,p) ^ in general, it admits as e — > 0+ a 
weak limit w(t) € Ai + (M. d , x K^), usually called Wigner measure (or semi-classical 
measure). The latter is known to give the possibility to describe in a "classical" 
manner the expectation values of physical observables via 

lim(ifj £ (t),Op £ (a)'i(j £ (t)) L 2 = // a(x,p)w(t,x,p)dx dp, 

£->0 I laid 
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where the physical observable Op e (a) is a self-adjoint operator obtained from the 
classical symbol a £ C£°(R!j5 x R^) through Weyl-quantization, see, e.g., p~3l[27] for 
a precise definition. 

Similarly to that, we were able to establish in |19| the existence of a limiting 
non-negative phase space measure /3(t) £ A4 + (R^ x R^), such that, after extracting 
an appropriate sub-sequence (denoted by the same symbol): 

/3 £ ^°> + /3 in L°°(R t ;M + (R d x x R^))weak*. 
If, in addition, ip E {t) is £- oscillatory, i.e. 

(1.10) sup {U £ (t)\\ L 2 + UeV^WIM < +co. 

0<£<1 

one can prove (cf. [191 Lemma 3.2]) that the limiting phase space measure /3(f) 
incorporates the classical limit of the position and current density in the sense that 

(1.11) P e {t,x) e ^ f /3(t,x,dp), 
and 

(1.12) J%t,x) £ ^\ + [ p(3(t,x,dp). 

Hereby the limits have to be understood in Ai + (M. d ,) weak*, uniformly on compact 
time-intervals J C If. Note that condition (jl.lOjl is satisfied for all t £ R, in view 
of (|1.2[) , the conservation of energy and our initial assumption (|1.3p (since for any 
V £ L°°(M. d ) we can be, assume without loss of generality, V(x) > 0). This is the 
reason to impose (|1.2p and (|1.3[) . throughout this work. In summary, the limiting 
Bohmian measure /3(f) yields the classical limit of the quantum mechanical position 
and current densities, by taking the zeroth and first moment with respect to p £ R d , 
analogous to the case of Wigner measures. 

While [B5] establishes several links between Wigner and Bohmian measures, it 
is not concerned with the dynamical properties of /3 e (f ) and its corresponding limit 
/3(f) (except for a short discussion on connections to WKB analysis before caustics, 
see [HI Proposition 6.1] and Remark 12.21 below). In contrast to that, the goal of 
the present work is do establish several dynamical result. The theorems obtained 
below can therefore be seen as independent from [19j , except for the basic existence 
and weak convergence results, recalled above. 

2. Main results 

The main objective of the current work is to analyze the dynamics of /3 £ (f). Note 
that the equivariance property strongly suggests that /3 e (f) satisfies the following 
nonlinear equation (of Vlasov-type) : 



d t ^+p- V X P £ - V x I V - 4^9^ ) • V P f3 e = 0, 



(2.1) 



p £ {t,x) = I f3 6 (t,x,dp), 
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subject to initial data (3 £ (x,p) given by (|1.6j) . Note that the system (|2.1[) can be 
written in one line as 

Due to the strong nonlinear nature of this equation and in particular due to possible 
singularities at points x £R d where p £ (t, x) = 0, it is a non-trivial task to rigorously 
prove that f3 e (t) furnishes a distributional solution to (|2.1j) . It will be one of the 
goals of this work to show that this is indeed the case. To this end, we shall derive 
appropriate bounds for all terms arising in the weak formulation of (|2.1j) after 
being evaluated at Bohmian measures. This rigorously establishes the existence of 
a nonlinear evolution equation for p E (t) in a similar spirit as the results of |12j for 
the quantum hydrodynamic system (jl.5p . More precisely, the first main result of 
this work is as follows: 

Theorem 2.1. Let V £ Cl(R d ;R) and % £ H 3 (R d ) with corresponding p e Q , J e 
defined by (|1.4[) . Then, for all e > 0, the Bohmian measure 

V P \t, x) 

is a weak solution of ([2T]) in V'(R t xl^xlj) and in V'([0, oo) x R d x R^) with 
initial data (11.61). 



In a second step we shall study the classical limit of equation (|2.1jl . To this end, 
let us recall the following definition used in [19]: 

Definition 2.1 (Mono- kinetic measure). A measure p £ M + (Rix Ep) is said to 
be mono-kinetic, if there exists & p £ Ai + (M. d ) and a function u(x) defined p - a.e., 
such that 

p{x,p) = p{x)S(p- u(x)). 

Mono-kinetic phase space measures define velocity distributions, which are given 
by the graph of u(x) and are thus particularly interesting for our purposes since, 
clearly, fi £ (t) is mono- kinetic by definition, its limit however, will not be in general. 

To get further insight, we pass to the limit e — > 0+ (after extraction of a sub- 
sequence) in the first three linear terms of equation (|2.1[) . This naturally leads 
to the following definition of a possible defect F(t, x,p) £ R d : Along a chosen 
sub-sequence {£„} ne N, let 

(2.2) T := ljm (div p (V x FJ/3 e )), in X>'(R t x R d x x R d ), 

such that 

A partial characterization of the defect J 7 , will be given in Section [5] In particular, 
it leads to the following result: 

Theorem 2.2. Let d = 1 and j3 e (t) solve (|2.1j) . In addition assume that at t = 0. 
the limiting measure satisfies 

fio(x,p) = w (x,p) = po(x)6(p - u (x)), 



(i 
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where wo(x,p) is the Wigner measure associated to iPq(x). Then, on any time- 
interval JCI| on which w(t) is mono-kinetic, it holds 

0(t,x,p) = w(t,x,p) = p(t,x)S(p-u(t,x)), 

in the sense of measures. 

In contrast to [T5], where we established certain sufficient criteria for ip e (t), 
yielding mono-kinetic Wigner and/or limiting Bohmian measures, this result di- 
rectly gives (3(t) = w(t) from the fact that the Wigner measure is mono-kinetic. 
It is well known, cf. jT2j [22] that the Wigner measure is generically mono-kinetic 
before caustics, i.e. before the appearance of the first shock in the (field driven) 
Burgers equation 



Note however, that there are situations in which w is of the form given in Theorem 
12.21 even after caustics, see e.g. Example 4 given in [12] and Example 1 in [27] , 
Both ot these examples furnish so-called point caustics, i.e. caustics where all the 
rays of geometric optics cross at one point. The result given above therefore shows 
that in some situations the limiting Bohmian measure j3{t) can indeed give the 
correct classical limit (for all physical observables) even after caustics. In addition, 
Theorem 12.21 generalizes p~9l Proposition 6.1], at least in d = 1 spatial dimensions. 

Remark 2.2. Equation (|2.3[) can be seen as the formal limit obtained from (|l.ip 
by means of WKB analysis. One thereby seeks solution to (jl.ll) in the following 
form 



with (real-valued) e-independent phase S and (possibly complex- valued) amplitude 
a 6 ~ ciq + ea\ + e 2 . . . , in the sense of asymptotic expansions. Plugging this ansatz 
into ([l.ip and neglecting terms of order C(e 2 ) yields (|2.3[) upon identifying u = V5, 
cf. [27] for more details (see also Section 6 of [19] where the connection between 
WKB analysis and Bohmian measures is briefly discussed). 

Finally, we shall consider the particular case where the initial data ipQ is a so- 
called semi-classical wave packet, see Theorem 12.31 below. The classical limit of 
Bohmian trajectories in this particular situation has been recently studied in [5J. 
There it has been shown that the Bohmian trajectories X s (t, x) converge (in some 
suitable topology) to X(t), the classical particle trajectory induced by the Hamil- 
tonian system 



One should note that the mathematical methods used in [9] are rather different 
from ours and that no convergence result for P E (t) is given, except for po = (and, 
as a variant, for a class of time-averaged Bohmian momenta). In comparison to 
that, the use of f3 e (t), together with classical Young measure theory, allows us to 
conclude the following result: 

Theorem 2.3. Let V G C^(R d ). In addition, assume that ipQ is of the following 
form 



(2.3) 



d t u + (u ■ V)u + VV(x) = 0. 



ip E (t,x) =a e (t,x)e tS{t ' x)/e 



(2.4) 




X = P, X(0) = x , 
P=-VV(X), P(Q)=p . 
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with given e -independent amplitude a G <S(R d ; C), such that |a(a;)| 2 > a.e. on M. d . 

(1) Then, the limiting Bohmian measure satisfies 

f3(t,x,p)=w(t,x,p) = \\a\\ 2 L2 S(x-X(t))S(p-P(t)), 

where X(t),P(t) are the classical trajectories defined by (|2.4I) . 

(2) Consider the following re-scaled Bohmian trajectories 

Y e (t,y) = X £ (t,x + V^y), Z £ {t,y) = P £ {t 7 x + ^y), 
where X £ (t,x),P e (t,x) solve (|1.8[) with initial data (xq,Pq) G R 2d . Then 

Y £ £ =^ X, Z £ ^ P, 

locally in measure on R t x R^ 7 i.e., for every 5 > and every Borel set 
Q C R t x R^ luzf/i finite Lebesgue measure TL, it holds 

(2.5) IimL({(t >V )Gfi: | p), Z%t, y)) - (X(t), P(t))\ > 6}) = 0. 

Let us compare now our results with the corresponding ones in [3]: Besides the 
fact that we are able to infer convergence of the Bohmian momentum P E (t), the 
topology we use is different from the one used there. Rephrased in our notation [5] 
proves that, for all T > and 7 > there exists some R < 00 such that, for all e 
small enough, 

(2.6) P pf) ({x G R 3 | max \X e {t, x ) - X(t)\ < R^}) > 1 - 7, 

where P p e is the probability measure with density = |?/>q(2;)| 2 . We see clearly that 
the comparison between the two results is not straightforward, as (|2.6|) measures the 
"good" points and ()2.5j) the "bad" ones. Moreover, (|2.6j) is more precise for finite 
times, whereas (|2.5[) doesn't require a-priori bounds on the time dependence and 
additionally implies almost everywhere convergence of sub-sequences, see Remark 
16.31 It would certainly be interesting to study the link between the two approaches 
more precisely. 

Remark 2.3. In view of Theorem l2.2l and Theorem [231 one might guess that w = j3 
only if both are mono-kinetic phase space distributions. This is consistent with the 
examples given in |19j but still wrong in general. Indeed, in the appendix of the 
present work, we shall construct a family of wave functions tp E for which w = (3 is 
absolutely continuous with respect to the Lebesgue measure on R^. 

The paper is now organized as follows: In the upcoming section we collect some 
basic mathematical estimates needed throughout this work. In Section [4] we estab- 
lish the fact that j3 e (t) indeed furnishes a distributional solution of (|2.1j) . In Section 
[5] the defect T will be partially characterized, yielding the proof of Theorem [221 In 
Section [51 the particular case of semi-classical wave packets will be studied includ- 
ing the proof of Theorem [231 Finally, in Appendix A we shall present and example 
in which w = /3 but not mono-kinetic. 

3. Static estimates 

In order to establish the weak formulation of (|2.1|) . we shall heavily rely on the 
following static, i.e. time-independent, estimate. 
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Proposition 3.1. Let ip £ C 2 (R d ), then it holds: 

|div (v|v| ® VM)| < d [ sup \d e d^d k 4>\ + i 



Im 



T 



sup \d e ipdjip\ , 



with j,t,k = 1,. . . ,d. 

Proof. We denote 9j := 9 X . and first compute 

vpdjip + tpdjip 



2(?/>vo 1/2 



which yields 



Using this wc can write 



(3-1) 



fc=i 



div (V|V| ® V|V|) = ^ Re cU fl^V + ^ ) . 



where each term in this series will be estimated separately (and in the same way). 
In order to handle terms in which the partial derivative dk acts on tp/tp we note 
that 



and we henceforth obtain 



dkip _ d k %pip _ ( dkTp _ <9fcV>\ ip 
ip ip 2 \ ip ^ / "0' 



i 



d e ipdjipd k — = -Im — — d e ipdjip— 







Using this on the r.h.s. of (|3.ip and summing up all terms yields the assertion of 
the lemma. 

□ 

In the upcoming analysis, we shall use the established estimate in the following 
form, where we denote by ||/||#i := ||V/||i2 the usual H 1 (R d ) semi-norm: 

Corollary 3.2. Fix e > 0. Then for any ip e £ H 2 oc (H d ) and for any test function 
if £ V(Ri x R d ), we have 



(3.2) 



je 

div (V^fff <g> V sff) <p[x,— 



dx < M £ < +oo, 



where p e , J £ are defined in (|1.4p . Explicitly, we find 



+ £rfi|V' e |!iji ( ml|V?/' e || i ji ( m sup / \ip{x,£)\dx, 

tm d JR d 

where fl C M~, denotes any compact set containing the support of >p(-,p) for all 
p £ R d . 

Proof. The proof follows directly from the estimate given in Proposition 13. II and a 
density argument in H 2 oc (M. d ). □ 
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In order to understand how M e behaves with respect to e we first note that due to 
the semi-classical scaling of the equation (jl.ljl ||V' E (*)l|jji is n °t uniformly bounded 
as e — > 0+ . In fact ip e (£) is e-oscillatory for all £ S R and thus each derivative scales 
like 1/e. Invoking the conservation laws of mass and energy, we have to use the 
re-scalcd semi-norms ||/|| ^ t := ||eV/||^2 in order to write M 6 in terms of uniformly 
bounded (w.r.t. e) expressions. Doing so, we find that M £ = 0(l/e 3 ) as e — > + . 

Remark 3.3. Note, however, that if tp e is of WKB type, i.e. ip e {x) = a(x)e lS( - x ^ £ 
with real-valued phase S(x) £ R and e-independent amplitude a, then the left hand 
side in (|3.2p obviously is bounded as e — > 0, whereas the right hand side blows up. 
In addition, it is easy to check that if one takes a superposition of two WKB states 
(such that VSi ^ VS2), the bound Q3.2p cannot be improved in general. 

4. BOHMIAN MEASURES AS DISTRIBUTIONAL SOLUTIONS OF (|2.1|) 

4.1. Mathematical preliminaries. Let us first note that the assumption V £ 
C^(R d ;R) is sufficient to ensure that, for each e > 0, the Hamiltonian operator 

H^ = -^A + V(x), 

is self-adjoint on D(H e ) = H 2 (M. d ) C L 2 (M. d ). It therefore generates a unitary 
(strongly continuous) group U £ {t) = e - ltHE / E on i 2 (R d ), which ensures the global 
existence of a unique solution ^ e (£) = U £ (t)ip of the Schrbdinger equation (jl.ll) . 
such that 

\\r(t,-)\\ L > = \m\L*- 

Moreover, since ipQ £ H 3 C D(H £ ), standard semi-group theory [53] implies ip e (t) £ 
D(H e ) = H 2 (R d ), for all £ e R. Clearly, this also yields that p e {t) £ L^R^R) as 
well as J £ (£) £ L 1 (M. d ;M. d ), for all £ G R and 

E £ (t) = E e (0) < +oo. 

In addition, since H e and U 6 (t) commute, and we obtain that 

H(i/T/V(£)|| L2(K£i) = \\(HT /2 M L * m . 

Since V £ L°°(R d ) the latter is equivalent to the n-th Sobolev norm 

ii/ik- :=ii(i + ier /2 )/n^, 

and we immediately infer the following result: 

Lemma 4.1. Under the assumptions of Theorem \2.1[ ip e (t) £ i£ 3 (R d ) for all (el. 

With this result in hand, we are sure to be able to apply the estimates established 
in Section [3] 

4.2. Weak formulation of (|2.1[) . In order to make sense of /? e (£) as a weak 
solution of (|2.1I) . the main problem is to understand the weak formulation of 
V x Vg • V p /3 e = div p (VgV ' x p e ). To this end, consider a class of compactly sup- 
ported test functions ip(t,x,p) = x(t,x)a(p) with X € C c °°(R t x R d ), a £ C c °°(R d ) 
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and compute 

POO P P 

V X V£ ■ V p (p(t, x, p)d(5 e (t, x, p)dt = 
x(t,x)V x V£(t,x) ■ V "a(u e '{t,x))p e \t,x) dx dt. 



JR d 



since, by definition, f3 6 (t,x,p) = p E (t,x)S(p — u e (t,x)) where denote u £ := '— 1 
i.e. the quantum mechanical velocity field. To this end, we recall that u e £ 
L 2 (R d , p e dx) . The following lemma then shows that this weak formulation indeed 
makes sense. 

Lemma 4.2. Let e > 0. For a £ C C °°(R^) and X G C c °°(E t x Rj*) we have 
\ X (x,t)\ \Va(u e (x,t))\ \V x VUx,t)\p e (t,x)dxdt < +oo. 



JWL d 

This result is the key for proving that Bohmian measures furnishes a distribu- 
tional solution of (12. ip . 

Proof. A simple computation shows that 

'Vp e ® Vp e 



P £ V^J = - VAp e - — div 

e 2 

= — V Ap e - e 2 div (Vy^ ® Vv^) . 
Inserting this into the weak formulation of V^Vg • V p /3 £ , we can estimate 

\x{x,t)\ \Wa(u £ (x,t))\ \V x V§(x, t)\p £ (x) dx dt 



e 2 r°° 



- Y J I \x(^t)\\^^ e (^t))\\^^P e \dxdt- 



■e 2 



\x(x,t)\ |Vcr(u £ (x,i))| | div (Vy/ff ® \dxdt 

JR d 



The two terms on the right hand side can the be treated as follows: By Lemma FTTl 
we have ip E £ H 3 (R d ) and thus VAp E is in L 1 (R), for all t £ M. Therefore 



\x{t,x)\ \Va(u e (x,t))\ \VAp £ \dxdt < +oo. 
On the other hand, inequality (|3.2[) directly yields (for any fixed e > 0) 

/ / |x(*.a:)| |Vcr(u 6 (a;,t))| |div(V\/p*® ^\TF) \dxdt < +oo, 
and the assertion is proved. □ 

As a final preliminary step, let us recall the classical notion of the push-forward 
for measures: Let po £ A4(M. d ) and / : M —> ]R d a measurable map. Then 
Mi = / $ Mo is called the push-forward of po under /, if for every cr G Co(M ), it 
holds: 

(4.1) / a(x)dp! = / o-(f(x))dp , 

JR d JR d 
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By using an approximation argument this condition can be relaxed in order to take 
into account test- functions a which are (only) integrable with respect to /Ui, but do 
not necessarily lie in Co, see Definition 1 and below in 

Remark 4.3. Indeed, the equivalent definition of a push-forward: 

Mi(0) = Mo(r 1 (0)), 

for all measurable sets C M. d and any measurable function /, allows one to choose 
indicator functions of measurable sets as test functions. Based on this observation, 
an approximation argument by simple functions leads to the fact that both sides 
of (|4.ip are finite and equal for all test functions a which are \i\ integrable [20] . 

The reason for using this slightly more general version of the usual push-forward 
formula will become clear in the proof given below. We first note that 



ip(t,x,p)df3 E (t,dx,dp)dt = / / ip{t,X E (t,x),u E {t,x))p E {x)dxdt, 

Jo Jm d 

by using the fact that P E {t, x,p) is the push-forward of the measure Pq(x)S(p~uq(x)) 
under the Bohmian phase space flow $f defined in (|1.7[) . 



Proof of Theorem POl Let ip g C£°([0, +oo) xRjx Rf) be a compactly supported 
test-function such that <p(t, x,p) = x{t, a; ) cr (p)- Then, the weak formulation of (|2.1[) 
reads 

OO f f 

o-{p){{d tX (t, x)+p- V xX {t, x))(i E {t, dx, dp)dt 

R 2d 

X (t, x)V x (V(x) + V§(t, x)) ■ V p a(p)l3 E (t, dx, dp)dt 

2d 

x(0,x)a(u E o (x))p E o (x) dx. 

First, consider the term involving V X V^: Having in mind the result of Lemma 14. 2\ 
we are allowed to consider x(t, x)W x V§(t, x) ■ Va(u E {t,x)) as a test-function which 
is integrable with repsect to p E (t, x). Thus, we can apply the push- forward formula 
P~T|) and infer 

f°° [[ x(t,x)V x V£(t,x)-V<T(p)[3 E (t,dx,dp)dt = 
o J Jm. 2d 

X {t, X E (t, x))V x V§(t, X s {t, x)) ■ Va(u E (t, X E (t, x)))p E {x) dx dt. 



In addition, the fact that J E = p E u E 6 L 1 (K d ) implies that the "test-function" 
a(u E (t, x))u E (t, x) ■ Vx(t,x) is integrable with respect to p E (note however, that 
u E (t,x) is not necessarily continuous). Thus we can again apply the (generalized) 
push-forward formula (|4.1j) to obtain 

a{p)p ■ V x x(t, x)/3 E (t, dx, dp) dt = 

a(u E (t,x))u E (t,x) ■ Vx(i, x)pq(x) dxdt. 
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All the other terms appearing in the weak formulation of (|2.1j) can then be treated 
analogously. Having in mind the ODE system (|1.8[) . we consequently arrive at 



((dtX + P- V x x)cr(p) - x^x{V + Vg) • V p a) P £ (t,dx,dp)dt 

^ ( X (t, X 6 (t, x))a(u e (t, x))) plidx) dt 
• d at 

X(0, x)<r(u e (x))p 6 (dx). 

This proves that the Bohmian measure fi e (t) furnishes a weak solution of (|2.ip in 
O'([0,oo) x M. d x x R^J) with initial data CCI]). The proof for V'(R t x R d x R*) is 
analogous. □ 

5. Study of possible defects 

Having established the fact that j3 £ is indeed a weak solution of (|2.1|) . we first 
rewrite the equation in the following form 

(5.1) d t p e + p ■ V x f3 e - V X V • V p /3 £ = div p (V x V§f3 £ ), 

where Vb is the Bohm potential defined in (jl.9p . Using the weak convergence results 
given in |19) we can pass to the limit on the left hand side of this equation (up to 
extraction of sub-sequences) in order to obtain 

(5.2) d t p + p ■ V x /3 — V x y • Vp/3 = T 

where the defect T is, as defined in (|2.2[) . In order to gain some information on T 
(and to prove Theorem I2.2[) . we first derive from (|5.2[) the following equation for 
the first moment of (3(t) with respect tope M. d : 

(5.3) d t J + div x p®pP{t,x,dp)-pV x V = -l pF(t,x,p)dp, 

where p, J denote the classical limits of p £ , J e given by (|l.lll) and (|1.12l) . Note that 
the integral on the ride hand side is well-defined since energy conservation implies 

\p\ 2 ^(t,x,dp) <C, 

where the constant C > is independent of e. This allows us to conclude that all 
terms on the left hand side of (|5.3|) are well-defined in the sense of distributions 
on R t x M.f and hence also the right hand side is. To proceed further we need the 
following result. 

Lemma 5.1. Let p e M+iR^ x R^) be such that 

\p\ 2 p(dx 7 dp) < +oo, 



p(x, dp) = p e M + (R d ), / pp{x, dp) = p(x)u(x), 

: M 

for some function u(x) £ R d defined p - a.e.. Then it holds 

p® p n(x, dp) > p(x) u{x) ® u{x), 
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with equality if and only if p(x, dp) = p(x)8(jp — u(x)). 

Proof. The proof follows directly from the fact that for p - a.e. x we have 

0< / {p — u(x)) ® (p — u(x))p,(x, dp) = / p®pp(x 1 dp) — pu®u. 



with equality if and only if p(x, dp) = p(x)S(p — u(x)). □ 
Using the result of this lemma we can rewrite (|5.3[) as 

(5.4) d t J + div x (pu®u) + pV x V = - / pJ 7 (t, x,p)dp - div x (pB), 

Jr* 

with a defect B(t, x) > 0. In addition we know that 

B{t, x) = 0, if and only if, f3{t, x,p) = p(t, x)5(p — u(t, x)). 

On the other hand, we can consider the equation for (3 E (t), take first the moment 
w.r.t. p and then pass to the limit e — > 0+: Multiplying (|5.1|) by p £ M. d and inte- 
grating yields the equation for the current density in the quantum hydrodynamical 
system p.5[) . i.e. 



(5.5) dt J £ + div ( ' F ® J£ ^ + p e VV = p £ WJ, 

where we have used the fact that 



p® p(3 e (t,x,dp) 



since f5 e is mono- kinetic by definition. Using this, we can define a defect C(t, x) > 



via 

(5.6) Urn ( J£0 E J ) = / P®Pp(t,x,dp)+C(t,x), 

where the limit has to be understood in V(R t xEJ). In addition, we know that 

e 2 _ 
p e WJ = — VAp £ - e 2 div {Vyfif® ^\fif) , 



2 

where the first term tends to zero as e — > 0+ by linearity. This consequently yields 
(5.7) d t J + div(pu ®u) + p\7V = - div(A + pB + C), 

where A{t, x) is defined by 



(5.8) e 2 (Vv^O VVp 6 ") -4 in V'{R t 



In summary, we have the following partial characterization of J 7 . 
Lemma 5.2. The defect J- defined in (|2.2[) satisfies 



(5.9) / pF{t,x,p)dp = -div(A(t,x) +C(t,x)), 

where A, C are given by (|5.8j) . (|5.6| . respectively. 
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In a last step, this can now be compared with the classical limit of the quantum 
hydrodynamic system (|1.5|) via Wigner measures. In [12] it has been shown that 

J £ (t,x) — ? J(t,x) := / w(t,x,dp) 



which satisfies 

(5.10) d t J + d\v{pu®u) + pVV = -divpT, 

with a temperature tensor T(t,x) > 0. The latter is found to be equal to zero, if 
and only if w(t, x,p) = p(t, x)5(p — u(t, x)). This can now be used as follows: 



Proof of Theorem \2.2\ Let d = 1. Having in mind that, by assumption, the initial 
limiting Bohmian and Wigner measures are equal /3q(x,p) = wo(x,p), the unique- 
ness of solutions, together with (|5.10|) and (|5.7|) . implies 

(5.11) pT = A + pB + C, in V(R t x R^), 

Since all terms on the right hand side are greater or equal to zero, we infer that 

uu(t, x,p) = p(i, x)S(p — u(t, x)), if and only if, A = B = C = 0. 

By definition, this implies that 

P £ WJ 0, 

as well as 

lim / p®p(5 £ {t,x,dp)= / p (g) p(3(t,x,dp) = pu <g> u. 



By Lemma f5.ll we conclude ft(t,x,p) = p(t,x)S(p — u(t,x)) and the assertion is 
proved. □ 

Remark 5.3. In dimensions d > 1 we can not conclude as before, since identity 
(|5.11|) has to be replaced by 

pT = A + P B + C + V, 
for some T> satisfying divl?(t,a;) = 0. 

6. Bohmian measures for semi-classical wave packets 

The theory of semi-classical wave packets is very well developed, see e.g. [T¥l [T51 
[211 [22] and the references given therein (see also [1] [6] for a recent application in 
the context of nonlinear Schrodingcr equations). It allows us to approximate the 
solution to (jl.lj) via 

(6.1) r(t 7 x) e ~° + u £ (t, x) = e~ d l\ (t, X ^ (t) ^j e WH-*(*))+S(*))A 

where P(t),X(t) solve the Hamiltonian system (|2.4p and S(t) is the associated 
classical action, i.e. 

S(t)= J*l\P( s )f-V(X(s))d s . 

The envelope function v(t,y) is thereby found to be a solution of the following 
e-independent Schrodinger equation (see also the proof of Theorem 12.31 below) : 

(6.2) id t v = -^A y v + -(Q(t)y,y)v, v(t = 0, y) = a(y), 
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where Q(t) := HessU(X(i)) denotes the Hessian of the potential V(x) evaluated at 
the classical trajectory X(t) and a € 5(M d ) is induced by the initial data i/jq given 
in Theorem 12.31 In other words, v(t, x) solves a linear Schrodingcr equation with 
time-dependent quadratic potential. Under suitable assumptions on V (satisfied by 
the hypothesis of Theorem 12. 3j) . one can show, see e.g. [TJ [14j [15l HI [21j [22], that 
u £ (t,x) approximates the exact solution ip £ (t,x) of in the following sense 

(6.3) \\r(t, •) - u £ (t, -)||x»(r*) < CiV^ 2 *, 

provided the initial data ip £ (0,x) is of the form given in Theorem 12.31 In other 
words, the approximation is valid up to times of order ln(l/e), at most. 

Remark 6.1. Note that in contrast to the WKB approximation, the coherent state 
ansatz does not suffer from the appearance of caustics (although it is sensitive to 
them through equation (|6.2j) where the caustics are somehow hidden) . In addition, 
it assumes that the amplitude concentrates on the scale yje. The latter has been 
shown to be a critical scaling in the theory of Bohmian measures, cf. the case 
studies in [19]. 

Proof of Theorem \2.3[ We first note that the solution to (|6.2j) satisfies \\v(t, -)\\l 2 = 
||a||^2 for all tel. Thus, the Wigncr transformation of u e (t,x) satisfies 

w £ [u £ ] w in L°°(mt;S'(Ri x M^)) weak*. 

The corresponding Wigner measure is well known, cf. jTTl [19] : 

w(t, x,p) = |M|| 2 S(x - X(t))S(p - P(t)). 

From the estimate (|6.3[) and the classical results given in [17] we conclude that the 
Wigner transformation of the exact solution w e [7/> e ] converges to the same limiting 
measure w, uniformly on compact time- intervals / C Rf 

In order to prove that w(t) = (3(t), we perform the unitary transformation: 
V> £ i— > v £ , defined via 

(6.4) r{t, x) = s- d ' A v £ U X ~^ (t) ) e^ p ^- x ^+ s ^l £ . 



Using this transformation, equation (jl.ip is easily found to be equivalent to 
1 

2" 



(6.5) id t v £ = --A y v £ + V £ (t,y)v £ , v £ (t = 0, x) = a{x) 



where V £ (t, y) is given by 

V £ {t 1 y) = - £ (V(X(t) + v/iy) - V(X(t)) - J!VV{X{t)) ■ y) . 

Obviously, for C 2 potentials V equation (|6.5[) converges to (|6.2p as e — > 0+. This 
together with sufficient a-priori bounds on v £ (t) yields the estimate (|6.3[) , cf. [5] for 
more details. On the other hand, using (|6.4[) . the Bohmian measure f3 £ (t) of the 
exact solution ijj e (t) can be seen to act on Lipschitz test-function ip e Cq(M.^ x R^) 
via 



(F(t)^) = I \v E (t,y)\\(x(t) + ^y,^l m (Y^l 

Jill \ V V (*! 



v) 



P{t) dy. 
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Using the Lipschitz continuity of tp we can estimate 



<p[X(t) + y/ey,y/eIm. 



< C^e \y 



Im 



v e {t,y) 

v £ (t,y) 



+ P(t))-cp(X(t),P(t)) 



for some positive constant C' v > 0. In view of this, we obtain 
W £ (t),v)- [ \v £ (t,y)\ 2 ^(X(t),P(t))dy 

\y\\v £ {t,y)\ 2 dy + ^ [ \v £ {t,y)\\Vv s {t,y)\dy 



< C v y/e 



<CWi|K(t)|| L2 (|| \y\v £ {t)\\ 



|Vv e (t)lU^ 



where the last inequality directly follows from Cauchy-Bunjakowski-Schwarz. In 
order to proceed further we need the following lemma. 

Lemma 6.2. Let V e Cj^(K d ). Then the solution of (|6.5[) satisfies 

\\\y\v B (t)\\ L2 <C U \\Vv £ (t)\\ V 2<c 2 , VieKTjci, 

where C\,Ci are some positive constants, independent of e. 

Proof of Lemma \6. 6 A In [5] it is shown that, if V is quadratically bounded, i.e. 
dW(x) e L°°, for all > 2, it holds: 



(6.6) 



\\\y\v £ (t)\\L* <Ci, \\\y\ 6 v £ (t)\\ L 2 <C 3 ,Vte 



It therefore only remains to show the estimate for Vv £ (£). This follows by consid- 
ering the energy corresponding to (|6.5[) . i.e. 



which satisfies 



E £ (t) = ^- [ \Vv £ (t,x)\ 2 dx+ [ V £ (t,x)\v £ (t,x)\ 2 dx, 

^-E £ (t)= f d t V £ (t,x)\v £ (t,x)\ 2 dx. 
at J M d 



Since 



\d t V £ (t,x)\ < \X(t)\\y\ 3 ^svp\d 3 V(X(t) + sVty)\, 
the assumption V € C 3 , together with (|6.6p . yields the desired bound on Vi/ E (i). □ 

Using the a-priori estimates established in Lemma 16.21 we obtain 



0. 



{P £ {t), v )- / \v £ {y,t)\^{X{t),P{t))dy 
In other words, we have that 

(6.7) P{t) = \\v(t, -)f L ^{x - X(t))S(p - P(t)), 

Having in mind that \\v(t, -)\\l 2 = IMIl 2 this proves assertion (1) of Theorem 
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In order to conclude assertion (2) of Theorem l2.3l we recall the following formula, 
stated in [HI Remark 3.8]. For all i G R and for all test-functions tp G C (R^ x R^), 
X G C (M t ) it holds 



X (t) // <p(x,p)P s {t,dx,dp)dt = 
: J JR 2d 

(6.8) = / X [t) [ ^(X £ (t 1 x),P £ (t,x))p £ (x)dxdt 

JR d 

x(t) ( p(x £ (t,x Q + V^y),P £ {t 7 x + V^y)My)\ 2 dydt, 



where in the second equality we set y = (x — xo)/y/s and recall that the initial 
density is given by 



p%{x) = £ 



- p -d/2 



X — Xq 

~7T 



Now, let the function ( e (t, y) be defined as 

C : Ri x R d y 3 (t, y) ^ (Y £ (t, y), Z 6 (t, y)) elfx R* 

where 

Y e (t, y) = X E (t, .T + y/Zy), Z £ (t, y) = P £ {t, x + Vh), 

are the rescaled Bohmian trajectories. The associated Young measure (or, parametrized 
measure): 

u t , v : t( x ij -) M + (Rt x Rp) ; (t, y) ^ u t , v (x,p), 
is defined, after extraction of a subsequence, via 

lim / / a(t,y,Q(t,y))dydt = // a(t,y,x,p)dw t , v (x,p) dy dt, 

J JS. t xtj J JWLt xM d J JR 2d 

for any test function a £ L 1 (R t x R d y ; C (R 2d )), cf. [H EH [26]. 

By passing to the limit e — > 0+ in (|6.8[) we find that (after the choice of an 
appropriate sub-sequence) : 

X(t) // tp(x,p)P(t,dx,dp)db= / x(t) / ¥>(x,p)ut, v (dx,dp)\a{y)\ 2 dy dt, 



and thus 

/3(t,x,p) = / \a{y)\ 2 uj t ,y{x,p)dy. 

Upon inserting (|6 . Tf) with \\v(t, -)\\l 2 = IMIl 2 > this implies 

Wt^Cx, jj) = u(t, y)S(p - P(t))S(x - X(t)). 
Since < v{y,t) < 1, |a(y)| 2 > 0, by assumption, and 

l«(y)l 2 rfy= / |a(v)lM*,tf)di/, 



for all t £ I, we conclude i/(t, y) = 1 a.e. with respect to |a(y)| 2 dy and hence 
w t ,„(a:,jj) = <J(p - P(t))6(x - *(*))• 

In other words, uit,y is found to be concentrated in a single point in phase space 
R^ x Rp. By a well known result of Young measure theory (see e.g. [T6l Proposition 
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1]), we know that the fact that &t,x is concentrated in a point is equivalent to the 
convergence of the re-scaled trajectories, i.e. 

locally in measure on R t x Rjj. □ 



Remark 6.3. In d = 1, assertion (1) of Theorem l2.3l directlv follows from Theorem 
12.21 since w(t) is obviously mono-kinetic. In addition, one should note that the 
established local in measure convergence implies (see e.g. [5J Section 13]) that 
there exists a sub-sequence {£Vi}neNi going to zero as n — > oo, such that 

yEn ™ x ^ z£n ™ p a e m ^ c M t X R d y . 

The rescaling of the position variable x = x + yfey in the Bohmian trajectories 
is critical for the coherent wave packet of Theorem 12.21 in the sense that it allows, 
at any time t £ R, to directly connect the Bohmian measure associated to the wave 
function with the Young measure of the Bohmian trajectories (with rescaled initial 
position). It seems that this cannot be done for any other rescaling of u £ (e.g. 
replacing yfe by £°). 

Acknowledgment. The authors want to thank the referees for helpful remarks 
in improving the paper and in particular for the short, direct argument given in 
the proof of Lemma 15.11 

Appendix A. An example with non mono-kinetic limiting Bohmian 

MEASURE 

In [19] we considered several different examples of ip £ and computed the corre- 
sponding limiting Bohmian measure /3 and the corresponding Wigner measure w. 
We found that in general w ^ (3 except in rather special situations. In fact, in all 
the examples given in |19j we find w = (3 only in the mono-kinetic case. Together 
with the results stated in Theorem l2.2l and Theorem 12 .31 this might yield the wrong 
impression that w and j3 can only coincide if they are both mono-kinetic phase 
space distributions. The following example will illustrate that this is in general not 
the case: 

Consider an e-dependent family of wave functions {u e }o<e<i given by 

u £ (x) = a e (x)e lS ^/ £ , 

where the amplitude a e reads 

a s {x)=e- d ' A p 1 ' 2 ( 



1/2 



with some e-independent profile p £ <S(R;R), i.e. smooth and rapidly decaying. In 
addition, let S £ Cb(R d ) H C^(M d /{0}), such that S{x) is even and has a cone-like 
singularity at x = 0: 



lim VS{Su) = x(w), £ § 



d-l 



for x £ C 00 ^- 1 ) 



Lemma A.l. Let u e be as given above, then 

] 

w 



f3(x,p) = w(x,p) = -rr I p(\y\)dy I 8(p - x(w))dw ® S(x) 
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Remark A. 2. To our knowledge this is the first example in which the projection 
of /3 = w onto position space is absolutely continuous with respect to the Lebesgue 



measure on lii d . 



Proof. We first note that, by assumption, 

|2 _ P ~d/2 n ( _M_\ ^0+ 



\ u *( x )\* =e -*l*p\JZj-j -^5(x) J p(\y\)dy. 

Moreover, it has been already computed in [12], that 

K™ + {™ e [u s ],<p) = Tg^ij J^p{\y\)dy ^(0, X (uj))duj, 

for any test function ip £ Co(K 2d ). Thus, we see that the Wigner measure uu(x,p) 
is of the form given above. 

It remains to explicitly compute the limiting Bohmian measure. To this end, we 
consider the action of (3 e onto any test function tp(x,p), i.e. 

(F,<p) = s- d / 2 J p (J^ ip{x, VS{x))dx 

) 

p(r)ip(e 1 / 2 rw 1 WS(e 1 / 2 rw))r d ~ 1 drdu;, 



is*- 1 Jo 

by setting y = rcu. It easily follows that as e — >• 0+: 

(P e ,cp)~f f p(r)r d - 1 ip(0,VS(e 1/2 rLu)drduj. 
Js^-i Jo 

Keeping r > 0, lu <G S d_1 fixed, we see that for e sufficiently small, 
By dominated convergence, we therefore conclude 

and the assertion is proved. □ 
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